The temperature influence on quantum tunneling in the spin-boson model 
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The existing studies of the temperature effect of the spin-boson model show that the tunneling 
splitting will increase with temperature, no matter how weak the couplings are between the bath 
and the spin. However, the small polaron theory said that in the weak coupling and low tempera- 
ture regime, the tunneling is, in fact, dominated by the diagonal transitions whereas this diagonal 
contribution to the tunneling should be faded away with increasing temperature. Taking advantage 
of the analysis originated from the Feynman- Vernon's influence functional theory, the influence on 
the tunneling by the phonon bath can be expressed as a product of the effects of every single phonon 
mode, which can be studied by numerical exact diagonalization. We find that, in the weak coupling 
and low temperature regime, all the spin-single-mode systems exhibit the same decreasing depen- 
dence of tunneling on increasing temperature. In turn, with the conclusion of influence functional 
theory, the tunneling of the spin-boson model should decrease and it is independent of the bath 
structure. In the strong coupling regime, however, the temperature effect reverses from suppressing 
to enhancing the tunneling with the increase of temperature. Discrepancies between the old theories 
and the small polaron theory are also explained. 

PACS numbers: 03.65.Yz, 03.65.Xp, 73.40.Gk 
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I. INTRODUCTION 

The spin-boson model which consists of a two-state 
system (TSS) coupled linearly to a phonon bath is a 
paradigm model of quantum dissipative systems. The 
study on this model has a long history and it still at- 
tracts a great deal of interest. 1-19 The basic question is 
how thermal fluctuations from the environment influence 
quantum tunneling of the TSS since all the real systems 
are at finite temperature. The answer to this question 
is important for the observation of macroscopic quantum 
phenomena and also for the realization of the quantum 
computer where the stability of a "qubit" against thermal 
fluctuations is critical. 20 Up to now, this issue has been 
addressed by various treatments^ ' 11 ! 12 ] 15 including nu- 
merical calculation, 7 and it seems that all got the same 
conclusion, i.e., quantum tunneling is enhanced by the 
temperature in low temperature regions. However, be- 
sides the different quantitative results, the preconditions 
to the conclusion are also different. Using a perturbation 
approach combined with a renormalization group (RG) 
treatment based on path integral, Bray and Moore firstly 
found that quantum tunneling increases with tempera- 
ture provided the coupling to the bath exceeds some crit- 
ical valued Numerical calculation based on path-integral 
also shown that the same conclusion can be found in the 
strong coupling regime but it was also stated that the sit- 
uation in weak coupling regime is qualitatively different 
from that in the strong coupling regimeZ On the other 
hand, the flow equation analysis^ and the variational 
calculation 9,15 suggested the enhancement of quantum 
tunneling happens in the weak coupling and the low tem- 
perature regions. Notably a quantitative result was de- 
rived and physics analysis on the enhancing mechanism 
was also provided in the variational calculation for the 
sub-Ohmic caseJ^ In fact, the quantitative result found 
in the sub-Ohmic case is in agreement with that found by 



Weiss^ such a T 1+s dependence was also considered as 
the universality in dissipative two-state system and was 
declared to be valid for all the coupling strength before 
the localization happens.— Nevertheless, physics analysis 
from the well-known small polaron theory tells a different 
story. 

It is known that the zero-biased spin-boson model is 
just a simple example of the polaron-phonon (or exciton- 
phonon) system, i.e., the two-site problem^ According 
to the small polaron theory^ the contribution to tunnel- 
ing has two parts at finite temperature, say, the diago- 
nal and non-diagonal contribution. In the weak coupling 
regime, the diagonal part makes the main contribution to 
quantum tunneling and hence quantum tunneling should 
decrease with increasing temperature in the low temper- 
ature regions, a result that is in confliction with some 
analysis of the spin-boson model mentioned above. The 
discrepancy between various treatments indicates that 
the consensus of temperature effect is still lacking and 
this is the main interest of the present work. 

In this paper, the temperature effect on quantum tun- 
neling in the spin-boson model is studied by analysis orig- 
inated from the influence functional theory developed by 
Feynman and Vernon. 4 One important conclusion of this 
theory is that the contribution of the whole bath can 
be expressed as a weighted integration over the contri- 
bution of single phonon modes; hence the temperature 
effect on quantum tunneling can be qualitatively known 
by studying the temperature effect in the single mode 
case which can be exactly done by numerical calculation. 
It is found that, in the weak coupling regime, all the 
phonon modes contribute the same temperature effect, 
i.e., quantum tunneling is suppressed by temperature, 
and the enhancing mechanisms by both the variational 
calculation 1 - and flow equation analysis^ do not exist in 
a real system. In the strong coupling regime, however, 
the temperature effect is reversed, in agreement with the 
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known results^ The organization for the rest of the pa- 
per is as follows. In the next section, the model and 
the explanation of our calculation method are presented. 
Numerical results for both T = and T ^ cases are 
given in Sec. III. Conclusion and discussion are presented 
in the last section. 



II. THE MODEL AND THE EXPLANATION 

The Hamiltonian of the spin-boson model is given by 
(setting h = 

h = --^-Vx +y^ y LQkbjh + ^y^gfe^i + h), (i) 

k k 

where — x,y,z) is the Pauli matrix, b k {b\) is 

the annihilation (creation) operator of the fcth phonon 
mode with energy to k and g k is the coupling parame- 
ter. The solution of this model depends on the so-called 
bath spectral function (density) which is defined as^ 
J(lu) — 7r J2k — For frequencies below a high- 

energy cutoff u) c , the spectral function has the power law 
form, 

J{u) = law'/wj- 1 , < w < lj c , (2) 

where a is a dimensionless coupling strength which char- 
acterizes the dissipation strength and the parameter s 
characterizes the property of the bath, i.e., < s < 1, 
s = 1, and s > 1 represent, respectively, the sub-Ohmic, 
Ohmic, and super-Ohmic cases. It should be noted that 
the solution of the model (1) is merely determined by 
the spectral function J(u>) comes from the fact that the 
bath degree of freedom can be integrated out as Gaus- 
sian integrals,^* 1 ^ an ideas originated from the influence 
functional theory.— 
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FIG. 1: Illustration of how the phonon bath is divided into 
a series of dynamically and statistically independent subsys- 
tems in continuum limit. Each subsystem represents a single- 
mode system for modes with frequencies lie in (u)k, Wk + duik]- 

As shown in Ref. |4j, the phonon bath can be considered 
as a distribution of single modes and the single modes can 
be treated as dynamically and statistically independent 



subsystems. An illustration of how the phonon bath with 
frequency lies in (0, uj c ] is divided in to a series of sub- 
systems, i.e., the single- mode systems, is shown in Fig. 1. 
According to the influence functional theory, 4 the total 
influence functional of the bath can be expressed as a 
product of that of each subsystem 

k 

where T k is the influence functional of the fcth mode. 
The validity of the above result is based on the assump- 
tion that the correlation between different modes can 
be ignored even when the interaction between the TSS 
and phonon modes is considered. This implies that the 
ground state of the whole system can be expressed as 

i $ °> = ^EwIMW)' (^(0)14(0)) = ***', 

(4) 

where |±) is the eigen-state of a z and |0?(O)) is the eigen- 
state of the fcth mode when it interacts with the TSS 
individually. Accordingly, one can find out the tunneling 
splitting 

A/A = <$o|<7 x |$o> = n rf fc(°)> ( 5 ) 

k 

where (4(0) = (0^(0)10^(0)) is the dressing factor con- 
tributed by the fcth mode. 

The above result can be generalized to the case 
by assuming that all the subsystems are still statistically 
independent when the interactions with the TSS are con- 
sidered. In the influence functional theory, a trace over 
each phonon mode is made in Eq. (3). 4 In the present 
case, the generalization can be done in two steps. Firstly, 
all the excited states of the whole system are suggested 
to have the same form as the ground state given in Eq. 
(4), that is, the nth excited state with eigenvalue E n can 
be expressed as 

|$ n ) = i=£|±>ni<^(n)>, (6) 

where (0fe(ra)|0i(n')) = Skk'^nn', then the tunneling 
splitting at T ^ can be found by 

At/Ao =^E„e-^«<$„|a x |$„) (7) 
= ^E„e-^-rL d k (n), 

where 

Z = 5>-*"S d k (n) = {^(n)\cf>^(n)), (8) 

n 

The assumption that each subsystem is still statistically 
independent means 

Z = l[z k , (9) 

k 
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where z k = ~^2 n e _,3e ™( fe ) is the partition function of the 
A:th mode interacting with the TSS individually. There- 
fore, one can find 

A T /Ao = TT D k (T), D k (T) = — V e- 0e ^d k (n). 
k Zk n 

(10) 

The above analysis shows that, by Eqs. (5) and (10), the 
contribution of the phonon bath to the tunneling splitting 
is a product of the contribution from every single mode. 
This is the starting point of our calculation. In the next 
section, the validity of such a simple product assumption 
(i.e., Eqs. (5) and (10)) will be verified numerically in a 
simple way. 

Although other analytical treatments, such as the vari- 
ational calculation and flow equation method, are not 
based on the influence functional theory, the expression 
of tunneling splitting found have exactly the same forms 
given in Eqs. (5) and (10). In fact, the dressing factor of 
a single mode can be expressed as the form of influence 
functional given in Rcf. 4, 

D k (T) = exp{-MT)}, (11) 

where ipk(T) can be considered as the dressing phase 
of the fcth mode. Since the solution should be deter- 
mined by the spectral function only, a natural choice of 
the dressing phase is 

MT) = glf(oJ k ,T), (12) 

which is just the form found by other analytical treat- 
ments, including the perturbation calculation, the varia- 
tional calculation, and the flow equation method in spite 
of their self-consistent forms 4 & 12 ' 15 i 16 This implies that 
the ideas of influence functional theory has been tacitly 
adopted by analytical treatments generally. It should 
be noted that, in the present work, the dressing factor 
Dk(T) is calculated numerically, and hence one can only 
get qualitative conclusion and no quantitative result can 
be found. 



III. TEMPERATURE EFFECT ON QUANTUM 
TUNNELING 

According to the above analysis, the influence of the 
bath on quantum tunneling can be known by solving the 
following single-mode system 

h = -^a x +iutfb + a z \{tf +b), (13) 

where X/u) is a dimensionless coupling parameter. The 
above Hamiltonian is solved by numerical diagonalization 
using the basis |±)<g){|n)}, here \n) is the eigenstate oib^b 
and n = 0, 1, 2, • • • . To find out the dressing factor D k (T) 
in Eq. (10), we shall take the first N eigenvalues {e n } 



and the corresponding eigenstates {<p n } for calculation 
in practice, here N is determined by 

P(e N ~ e x ) > L, (14) 

and L is a fixed number to control the calculation error. 
It is found that the result becomes iV-independent when 
L > 20. Since the sign of (f n \S' x \(p n ) is not relevant for 
the calculation of the total tunneling splitting, we have 

1 N 

£> fe (T)~-^ e-^\(cp n \a x \cp n )\, (15) 

Z n=l 

where z ~ e _,9e ™ . We have compared numerical 

result with other approach using a different basis,— and 
the result is the same within the calculation error. A de- 
tail discussion on numerical diagonalization of the single 
mode Hamiltonian can be found in Ref. 

In the following, we shall present a simple verification 
of Eqs. (5) and (10) by numerical calculation. This is 
done in the following way. Firstly, the numerical solu- 
tions for two given modes with frequencies lo\ and w 2 and 
coupling parameters Xi/wi and A 2 /w 2 are found separ- 
atively, and the tunneling splitting of the ground state 
Ai i2 and the temperature dependence of Ai j2 (T) can be 
found by Eq. (15). Secondly, the numerical solution for a 
two-mode system with the same frequencies and coupling 
parameters are done in the same way, and the tunneling 
splitting of the ground state Ai + 2 and the temperature 
dependence A 1+2 (T) can be found. According to Eqs. 
(5) and (10), one should have 

A 1+2 = A 1 A 2 , and A 1+2 {T) = A 1 (T)A 2 (T). (16) 

Typical result for the ground state (i.e., the case of 
T = 0) is shown in Fig. 2. It is seen that the valid- 
ity of the above equation depends on parameters uj\/uj2, 
wi, 2 /Ao, and X/uii^- However, one can find that the 
above equation holds quite well within the calculation 
errors when X/uJ\ t 2 < 0.1. The result for the case of 
T 7^ is shown in Fig. 3. One can see that, in this 
case, the validity of the above equation also depends on 
the temperature region. The Ai +2 (T) = Ai(T)A 2 (T) 
holds for X/u)\ t 2 < 0.1 in the temperature regions with 
ksT ' I Ao < 0.1. Since in a real thermodynamical sys- 
tem, the coupling parameter of each single mode should 
scale with 1/N and N — + 00, one can conclude that Eqs. 
(5) and (10) hold in the low temperature region with 
k B T/A < 0.1. 

The above analysis shows that numerical solutions to 
the single-mode systems can serve as the solid basis for 
understanding the effect of the whole bath. As we have 
shown, tunneling splitting found by analytical treatments 
with various approximations have the form given in Eq. 
(11), this implies that numerical solution of the single- 
mode system can be used as a touchstone to test qualita- 
tively the analytical results. Let us first see the result of 
the ground state. Figure 4 shows the dependence of the 
tunneling splitting of the ground state on the coupling 
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FIG. 2: Numerical ver- 
ification of the simple 
product assumption in 
the case of T — and 
Ao/lji = 1/3. The 
deviation between A1+2 
and A1A2 appears as 
the coupling strength 
increases but both A1+2 
and A1A2 show sim- 
ilar A/wi^-dependence. 
In the weak coupling 
regime with A/ui,2 < 
0.1, it is found that 
Ai+2 = A1A2 holds 
well within calculation 



parameter for various frequencies. For comparison, we 
also present the result by adiabatic approximation, i.e., 
A a /A = (0±|0 T ), where |0±> = exp{=F(A/w)(St -6)}|0> 
is the displaced-oscillator ground state. It shows that the 
numerical result agrees with the adiabatic approximation 
for Ao/w < 0.1, but discrepancy becomes obvious when 
Ao/w > 1/3 and the lower the frequency, the larger the 
deviation. The result shown in Fig. 4 provides a quantita- 
tive explanation for the fact that the usual Debye- Waller 
reduction by adiabatic approximation on tunneling often 
appears to be by too target As one can see from Fig. 4, 
in the weak coupling regime with X/lu < 0.1, the low fre- 
quency modes with frequency Ao/a-> > 5 have almost 
no effect on the tunneling splitting, while the Debye- 
Waller factor tells that the low-frequency modes lead to 
localization ^ 15 ' 17 As a matter of fact, one goal for ana- 
lytical treatment is to find an appropriate dressing factor 
better than the usual Debye- Waller factor. Considering 
the dressing phase given in Eq. (11), the result by the 
flow equation analysis at T = isi2 



MT = 0) = -2- 



9l 



A 2 ' 



(17) 



where A is the tunneling splitting of the whole system 
and was denoted as Aqo in Ref. Il2j. Interestingly, this re- 
sult is in agreement with the perturbation calculation up 
to the first order approximation with A is replaced with 
Ao<£ According to the above result, the bath is demar- 
cated by a frequency oj* — A(or lj* = Ao in perturbation 
calculation), the modes with frequency larger than ui* 



will suppress tunneling splitting, while the modes with 
frequency lower than to* will enhance tunneling splitting 
since the dressing phase changes its sign. However, our 
numerical calculation find no evidence for the existence 
of uo* . We find that all the phonon modes lowers tunnel- 
ing splitting in the coupling regime < (X/ui) < 2.0 at 
T = 0. On the other hand, the variational calculation 



gives^ i 15 ' 19 



MT = o) 



9l 



K + A) 2 



(18) 



Comparing with the usual Debye- Waller factor, i.e., 
tpk{T = 0) = — 2g\juj\, the above result effectively 
shifts the frequency origin from lu = to lj = A, viz, 
the contribution of the low frequency modes with fre- 
quency uj < A is ignored, an ideas similar to adiabatic 
renormalizationfM As one can see in Fig. 4, such a treat- 
ment is reasonable since the low frequency modes have 
less effect on tunneling splitting in the weak coupling 
regime. This explains the fact that the variational ground 
state is a good approximation to the true ground state 
in the weak coupling regimes. 1 ^ However, the variational 
calculation fails to trace the phonon-induced localization 
in the strong coupling regimeJ^ 

Now we turn to see the temperature effect on tunnel- 
ing splitting. Temperature dependence on the dressing 
factor for the single-mode system is calculated by Eq. 
(15) and some typical results are shown in Fig. 5. In the 
weak coupling regime, we have dAk(T)/dT < for all 
the phonon modes. For a given coupling strength with 
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FIG. 3: Numerical ver- 
ification of the simple 
product assumption 
in the case of T / 
with A/a>i,2 = 0.1 and 
Ao/cji = 1/3. The devi- 
ation between Ai+2(T) 
and Ai(T)A 2 (T) 
becomes obvious as 
temperature increases 
but both Ai +2 (T) and 
Ai(T)A 2 (T) always 
show similar tempera- 
ture dependence. In the 
low temperature regions 
of fcsT/Ao < 0.1, it is 
found that Ai +2 (T) = 
Ai(T)A 2 (T) holds 
well in the weak cou- 
pling regime with 
A/wi,2 < 0.1. 
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FIG. 4: Dependence of the tunneling splitting of the ground 
state on the coupling parameter (A/a;) for various frequen- 
cies. For comparison, the result by adiabatic approximation 
is shown in solid line. The numerical result agrees with the 
result by adiabatic approximation only when Ao/cj < 0.1, 
and the lower the frequency the larger the discrepancy. 
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FIG. 5: Temperature dependence on the dressing factor for 
the single-mode system in the case of (A/w) = 0.1 for various 
frequencies. We have dAk(T)/dT < for all the phonon 
modes. In the low temperature region with /cbT/Ao < 0.1, 
the phonon modes with frequency Ao/u> < 0.1 show almost 
no observable temperature effect. 



(A/oj) < 0.1, the lower the frequency (i.e., the larger 
Ao/w), the lower temperature where tunneling splitting 
starts to decrease. This is easy to understand since the 
contribution of the excited states appears at a higher 
temperature for a larger frequency. In the low tempera- 



ture region with ksT/Ao < 0.1, the phonon modes with 
frequency Aq/uj < 0.1 show almost no observable tem- 
perature effect. By Eq. (10), this result shows that, in 
the weak coupling regimes, tunneling splitting decreases 
with increasing temperature in the low temperature re- 
gion with fcsT/Ao < 0.1. It should be noted that, since 



6 



all the modes have the same temperature dependence, 
the result of whole bath according to Eq. (10) is qual- 
itatively independent of the weight of different modes, 
namely, it holds for the sub-Ohmic, Ohmic, and super- 
Ohmic cases. Obviously, this conclusion is in confliction 
with some known result a 2 ' 9 ' 11 ^ 12 ' 15 and explanations will 
be given below. 

According to the variational calculation^^ tempera- 
ture enhancing mechanism comes from the result that 
the dressing factors for frequency modes with frequency 
ksT < lo < A decrease with increasing temperature. 
We have searched for these enhancing modes in numeri- 
cal calculation, however, no evidence for the existence of 
such enhancing modes can be found in the weak coupling 
regime. The searching work is done in the following way. 
We perform a variational calculation for a single-mode 
system given in Eq. (13) by the same way as done in 
Ref. [HI then a self-consistent equation for the tunneling 
splitting of the single-mode system can be found 

A„(T)/A - exp{-(A/co) 2 / 1 (T, w )}, (19) 

where 

f(Tuj ) = coth(/V2) 

/U ' ' {u> + A„(T)tanh[/3A^(T)/2]coth(/3w/2)} 2 ' 

(20) 

from which the temperature dependence on tunneling 
splitting can be obtained numerically. It can be verified 
that, in the weak coupling regime (A/w) < 0.1, the vari- 
ational ground state is a good approximation to the true 
ground state and the temperature-enhancing modes do 
exist as predicted in multi-mode case,— that is, we have 
dA v (T)/dT > when k B T < lo < A v . Some typical re- 
sults are shown in Fig. 6, where one can see that, for the 
high frequency modes with A /lo < 1, the temperature 
dependence is qualitatively the same for both the vari- 
ational calculation and numerical calculation that tun- 
neling splitting decreases with increasing temperature. 
However, in the low frequency region with Aq /lo > 1, 
the result by variational calculation shows opposite tem- 
perature dependence to numerical result. This result in- 
dicates that the enhancing modes could be an artifact 
of the variational calculation in the case of T ^ 0. It 
should also be noted that another related prediction, the 
re-entrance phenomena^ has been shown to be impos- 
sible in a real system^ 

In the flow equation analysis, the enhancement of tun- 
neling splitting by temperature comes from the existence 
of a frequency boundary lo* = A. As one can see from 
Eq. (17), those modes with frequency lo < lo* = A al- 
ways enhance the tunneling splitting, while in the low 
temperature region with ksT < Ao, the main contri- 
bution of the whole bath comes from the modes with 
frequency < lo < A, this leads to the enhancement 
of tunneling splitting by temperature as shown in Ref. 
Il2l . However, as we have shown, there is no evidence for 
the existence of w*, hence the enhancing mechanism in 
this case does not exist. We believe that, in the case 



of T = 0, the enhancement of the low frequency modes 
by Eq. (17) can compensate for the over- estimation of 
the usual Debye- Waller reduction of the high frequency 
modes, which leads to a good final result. In the case of 
T^O, however, the fake enhancing mechanism can lead 
to a wrong temperature dependence in the low tempera- 
ture regions. 

In the above analysis, temperature enhancing mecha- 
nisms in the weak coupling regime by both the variational 
calculation and flow equation analysis have been ruled 
out. Now we turn to see what happens in the strong 
coupling regime. In this case, it is worth noting that 
the analysis in the strong coupling regime may not be 
as precisely as that in the weak coupling regime. This 
is because, in the strong coupling regime, the assump- 
tion of simple product (e.g., Eqs. (5) and (10)) does not 
hold exactly. However, our numerical result shows that, 
in the whole coupling range we studied, Ai + 2(T) always 
shows similar temperature dependence as Ai(T)A2(T), 
and furthermore, to the first order approximation (i.e., 
by omitting the correlation between different modes), the 
effect of the whole bath can be considered as a weighted- 
integration over different modes, hence the temperature 
effect in the strong coupling regime can be approximately 
traced in the same way as that in the weak coupling 
regime. In the following, it will be shown that the 
temperature enhancing mechanism appears in the strong 
coupling regime. 

According to Eq. (15), if the tunneling splittings of 
the first several excited states are larger than that of the 
ground state, the resulted Afc(T) will increase with tem- 
perature, otherwise A& (T) will decreases with increasing 
temperature. In other words, one can see the tempera- 
ture dependence by examining the spectrum of tunneling 
splitting. In the weak coupling regime and for all the fre- 
quency ranges we studied, K^nlcM^n)! decreases with 
increasing n for n < 5, leading to a decreasing tunnel- 
ing splitting with temperature; however, as the coupling 
parameter (A/cj) increases to some critical value, which 
is frequency-dependent, K^rilcM^™)! turns to increase 
with n and, of course, the temperature dependence re- 
verses. Typical results are shown in Fig. 7, where one can 
see the variation of the spectrum of tunneling splitting 
and the reversed temperature dependence as the coupling 
strength exceeds some critical value. Since the reversed 
temperature dependence in the strong coupling regime is 
found for all the frequency range we studied, one can con- 
clude that tunneling splitting increases with temperature 
in the spin-boson model provided that the coupling to the 
bath is strong enough, confirming the results having been 
found before££ Such an evolution of temperature de- 
pendence with coupling to the bath is in agreement with 
physics analysis of the small polaron theory. According 
to the small polaron theory^ the contribution to the 
tunneling has two parts: the diagonal and non-diagonal 
contribution. The diagonal contribution decreases with 
increasing temperature and can only survive in the weak 
coupling regime. On the other hand, the non-diagonal 
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FIG. 6: Temperature 
dependence of tunneling 
splitting in single-mode 
case by the variational 
calculation (shown as 
dot-line) and numerical 
calculation (shown 
as solid line). For 
high frequency with 
Aq/u) < 1, temperature 
dependence is qualita- 
tively the same for both 
cases that tunneling 
splitting decreases with 
increasing temperature. 
However, in the low 
frequency region with 
Ao/lj > 1, temperature 
dependence becomes 
opposite for both 
cases, showing that 
the enhancing mode 
is an artifact of the 
variational calculation. 



contribution increases with the temperature and only be- 
comes important when the coupling strength is not too 
low. In the weak coupling regime and low temperature 
region, the main contribution to the tunneling is the diag- 
onal part, which decreases with increasing temperature. 
In the strong coupling regime, the diagonal contribu- 
tion is suppressed and non-diagonal part makes the main 
contribution to the tunneling, hence tunneling splitting 
will increase with temperature. In the medium coupling 
range, a coherent-incoherent transition can happen when 
temperature increases as shown in Ref. l24l where one can 
also find a detail analytical analysis on temperature de- 
pendence of tunneling splitting for single-mode system 
by using the small polaron theory. The result shown in 
Fig. 7 provides an alternative explanation for the results 
found in Refs. [H and 0- 

IV. CONCLUSION AND DISCUSSION 

In the previous sections, the temperature effect on tun- 
neling splitting of the TSS is studied by a theory origi- 
nated from the influence functional theory. It is shown 
that the effect of the whole bath can be expressed as 
a simple product of every single mode, then the single- 
mode system is solved by numerical diagonalization, from 
which the effect of the whole bath can be found. In 
the weak coupling regime, all the phonon modes have 
the same effect on tunneling splitting of the TSS, i.e., 
coupling to any phonon modes suppresses the tunneling 
splitting, and the suppressing effect increases with tem- 



perature. The result indicates that tunneling splitting 
will decrease with increasing temperature for the sub- 
Ohmic, Ohmic, and super-Ohmic cases. Temperature 
enhancing mechanisms by both variational calculation 
and the flow equation analysis have been ruled out by 
numerical analysis. The reversed temperature effect on 
tunneling splitting in the strong coupling regime is also 
explained. Our analysis shows that the temperature ef- 
fect in the spin-boson model can be understood based 
on the small polaron theory. The present work can help 
to resolve the discrepancy between various treatments. 
However, in the present stage, it is not totally clear about 
the conclusion found in Refs. Q and [ll|. Here we just want 
to emphasize two points. First of all, the temperature ef- 
fect on quantum tunneling in the spin-boson model is 
hard to deal with since the usual Debye- Waller factor 
tends to zero even at T = for both the Ohmic and 
sub-Ohmic case. This implies that, the diagonal contri- 
bution to tunneling dies out by the small polaron theory 
and hence calculations done along this line is impossi- 
ble to include the diagonal contribution. As we known, 
temperature dependence on the diagonal contribution is 
opposite to that of the non-diagonal contribution. Ac- 
cordingly, it is possible to get a correct temperature de- 
pendence only when the calculation is made by a better 
approximation than the small polaron theory. Unfortu- 
nately most analytical calculations are done in similar 
way to the small polaron theory. By comparing the main 
result (see Eq. (9) in Ref. GH and EqJ_3.36) in Ref.Q]), it 
is seen that the result found in Ref. [ll| is similar to that 
found along the line of small polaron theory. 1 This seems 
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FIG. 7: Spectrum of tun- 
neling splitting for two typ- 
ical frequencies of Ao/cj = 
1/3 (up) and Ao/lu = 3 
(medium) in the weak (left) 
and strong (right) coupling 
regime(the line is guidence for 
the eyes). In weak coupling 
regime, (tp n \a x \ip n ) | decreases 
with n, while in the strong cou- 
pling regime, the dependence 
reverses. The resulted temper- 
ature dependence on tunneling 
splitting is shown in the lowest 
part, where one can see the tem- 
perature dependence reverses in 
the strong coupling regime. 
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to suggest that the diagonal contribution might be omit- 
ted in the calculation of Ref. [HI. The second point is the 
difference between the TSS system and the double-well 
system. In a double- well system, there are excited states 
inside each well and the tunneling splitting increases ex- 
ponentially with the excited state . 25 ' 26 The temperature 
dependence on tunneling splitting can be known without 
knowing the bath effect on tunneling splitting of each 
level. As the temperature increases, the occupation prob- 
ability of the excited states increases, which makes the 
tunneling splitting increases with temperature ! 25 ' 26 On 
the other hand, for a spin 1/2 system, there is not excited 



state for each direction. The temperature dependence on 
quantum tunneling depends on the spectrum of tunneling 
splitting when coupling to the bath. Unfortunately, even 
the bath effect on quantum tunneling in ground state is 
not exactly known in the spin-boson model. Some dis- 
cussion on temperature effect given in Ref. is for the 
double-well system, but not for the spin-boson model. 
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